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Abstract. Extremization of the Boltzmann-Gibbs (BG) entropy Sbg = —k J dx p(x)lnp(x) under ap- 

propriate norm and width constraints yields the Gaussian distribution pg(x) oc e . Also, the basic 
' solutions of the standard Fokker-Planck (FP) equation (related to the Langevin equation with additive 

noise), as well as the Central Limit Theorem attractors, are Gaussians. The simplest stochastic model 
with such features is N — > oo independent binary random variables, as first proved by de Moivre and 
Laplace. What happens for strongly correlated random variables? Such correlations are often present 
in physical situations as e.g. systems with long range interactions or memory. Frequently (/-Gaussians, 
p q (x) oc [1 — (1 — q)l3x 2 ] 1 ^ 1 ~ q ^ [pi(x) — pg(x)\ become observed. This is typically so if the Langevin equa- 
tion includes multiplicative noise, or the FP equation to be nonlinear. Scale-invariance, i.e. exchangeable 
binary stochastic processes, allow a systematical analysis of the relation between correlations and non- 
Gaussian distributions. In particular, a generalized stochastic model yielding g-Gaussians for all (q 7^ 1) 
was missing. This is achieved here by using the Laplace-de Finetti representation theorem, which embod- 
ies strict scale-invariance of interchangeable random variables. We demonstrate that strict scale invariance 
^ • together with g-Gaussianity mandates the associated extensive entropy to be BG. 

OO 

CO ■ PACS. 05.20.-y Classical statistical mechanics - 02.50.Cw Probability theory - 05.90.+m Other topics 

in statistical physics, thermodynamics, and nonlin. dyn. systems - 05.70.-a Thermodynamics 
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One of the cornerstones of statistical mechanics is the in thermal equilibrium with a thermostat i.e., if the inter- 
CTN . functional connection of the thermodynamic entropy with actions between its elements are sufficiently short-ranged, 
■ the set of probabilities {p{\ of microscopic configurations, or inexistent. Furthermore, this important probabilistic 
For the celebrated Boltzmann-Gibbs (BG) theory, this form (ii) maximizes the (continuous version of the) en- 
functional is given by tropy Sbg — — J dx p(x)hi[p(x)] under the basic con- 
straints of normalizability and finite width; (ii) consti- 
tutes the exact solution of the simplest form of the (lin- 
Sbg = —k 2_^Pi\iiPi , (1) ear) Fokker-Planck equation (i.e. corresponding Langevin 
1=1 equations with purely additive noise); (iii) is the TV — > 00 
where N is the total number of microscopic states which attractor of the appropriately centered and scaled sum 
are compatible with the information available about the of N ^dependent (or weakly correlated m an appropn- 
system. This powerful connection is in principle applicable ate se f e ) discrete ° r continuous random variables whose 
to a vast class of relevant systems, including (classical) dy- seco * d moment is finite (Central Limit Theorem). The 
• , 1 -it 4. • simplest probabilistic model which realizes these paradig- 
namical ones whose maximal Lyapunov exponent is posi- v 1 it 

five warranting strong chaos, hence mixing in phase space, matic features is a set of TV independent equa binary ran- 

hence ergodicity (Boltzmann, in some sense, embodied all dom ™bles (each of them taking say the values and 1 

these features in his molecular chaos hypothesis). Within wlth Probability 1/2). The probability of having, for fixed 

this theory exponential distributions, i.e. exponential sta- N > n Vs is S lven b ^ m (N- n y. 2 ■ Its limiting distribution 

tistical factors, emerge naturally. is, after centering and scaling, a Gaussian (as first proved 

In particular, the Gaussian distribution p G oc e~^ b y de Moivre and Laplace), and its (extensive) entropy is 

with P > is found to be the (i) velocity (Maxwell) dis- the one of BG -, since Sbg(N) = Nkln2. 
tribution of any classical many-body Hamiltonian system 
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What happens with the above properties when these 
binary random variables are not independent and the cor- 
relations between those variables (homogeneously spread 
over the system) are strong enough? There is in principle 
no reason for expecting the relevant limiting distribution 
to be a Gaussian, and the entropy which is extensive (i.e., 
S(N) cx N for N > 1) to be S B g- The purpose of this 
paper is to answer such questions for a mathematically in- 
triguing class of correlated processes - which can be inter- 
preted as prototypical mean field models - and therefore 
are relevant in natural, artificial and even social systems. 
Let us discuss the frequently occurring q-Gaussians, a nat- 
ural generalization of the Gaussians, defined as 



(2) 



O 7 

p q (x)<x[l-(l-q)px 2 ]W-ti \p 1 (x)=p G (x)} , 

where x £ R for q < 3 (for q > 3, normalizability is 
lost), and x 2 < 1/[(1 — q)(3] for q < 1. (i) q-Gaussians 
appear, e.g., as the exact solutions of paradigmatic non- 
Markovian Langevin processes [TJ[2] which lead to inhomo- 
geneous linear [3] , or homogeneous nonlinear [4, 5] Fokker- 
Planck equations, (ii) q-CLT attractors are q-Gaussians 
[6]. (iii) The extremization of S q with norm and finite 
width constraints yields q-Gaussians. S, is a (generi- 
cally nonadditive [8]) generalization of BG entropy, namely 



Sq — k 



1 - / dx \p{x)} q 



1 



(q G R; Si = S BG ). (3) 



Usually the extensive entropy of systems is given by the 
Boltzmann Gibbs entropy Sbg which is identical with S q 
for q = 1. It is noteworthy at this point that there exist 
systems with distribution functions (which are in general 
not q-Gaussian) depending on the system size N in such a 
way that S bg is not the extensive entropy any more while 
S q becomes the extensive entropy for specific values of 
q < 1 denoted by q cnt , i.e., S qmt (N) cx N (N > 1), [HU]. 
As is well known, for all standard short-range-interacting 
many-body Hamiltonian systems, we have q ent = 1, i.e. 
Sbg, which is identical to Si, is extensive, (iv) Numerical 
indications [llj for the distributions of velocities in qua- 
sistationary states of long-range Hamiltonians [12] suggest 
g-Gaussians. Different interpretations of the situation are 
given in [131114) . (v) Further, experimental and observa- 
tional evidence for g-Gaussians exists for the motion of bi- 
ological cells [TSlfTrj] , defect turbulence [T7] , solar wind [TBI 
[L9] . cold atoms in dissipative optical lattices [20], dusty 
plasma [21], among others. Numerical indications are also 
available at the edge of chaos of unimodal maps [55]. In 
[24: a specific model for correlated binary random variable 
was shown to converge to q-Gaussian distributions in the 
thermodynamic limit for q < 1. Yet, no mechanism was 
given of how such models can be generated from general 
principles. This will be done in this paper. In particular, 
it will be shown how the model given in [24] can be gener- 
ated and how by the same means q-Gaussian models can 
be generated for 1 < q < 3. 

In the following we consider binary exchangeable sto- 
chastic processes, with correlated elements say from x G 
{0,1}. Exchangeable means that the Appoint probabili- 
ties 



Pn(x\, X2, ■ ■ ■ , xjsr) of the process are totally symmetric 
in its arguments for all N and that p^ can always be 
obtained by marginalization of pn+i- In particular, the 
probability of a specific sequence (xi,X2, ■ ■ ■ ,xn), a rni- 
crostate, does not depend on the order of binary events, 
but only on the number n of events in the state Xi — and 
N — n events in the state Xj = 1. Following the notation in 
[24] we denote this probability as r% • There are ( ) such 
micro-states (where the order of the Xi is exchanged) , and 
the probability of finding any situation with n events in 
one state is given by 



P% = 



N' 



n 



.N 



N 
n=0 



(4) 



JV-l 



Total symmetry of the micro-states, VW, implies r. 
T"n + I'n+i' sometimes referred to as the Leibniz trian- 
gle rule [9], or scale-invariance of the distribution. Scale- 
invariance means here that at every 'scale' N the rela- 
tion p N (xi,x 2 , ■ ■ ■ ,x N ) = J2x N+ i Pn+i(xi,x 2 , ■ • ■ ,x N+ i) 
holds (where the sum runs over the possible discrete val- 
ues of the random variable), i.e. the distribution at a lower 
scale can always be obtained by marginalization of higher 
scales. What does this mean for physical systems? If we 
look at the binary stochastic process in terms of coin 
tosses, this means that the transition probabilities that 
bias the iV'th coin toss depend only on the number of 
tails and heads thrown after N — 1 tosses and not on the 
particular history of coin tosses. If the coin tosses steer a 
random walk this means that the next step of the random 
walk is biased by the elapsed time (number of coin tosses 
N) and the distance the random walk has covered to its 
origin (N — 2n) but not on the particular path the ran- 
dom walk has taken. In all Hamiltonian systems where the 
interactions are sufficiently long-range and strong enough 
for mean field solutions to be exact or close to exact to- 
tally symmetric correlations become physically relevant. 
Clearly, when mean field solutions become exact correla- 
tions between the state variables become totally symmet- 
ric. For binary state variables, i.e. Ising like spin systems, 
this implies that the conditional probability of one arbi- 
trary spin being in one particular state conditioned on 
the state of all other spins only depends on the number of 
spins in one or the other state and not on the particular 
spin configuration of the system. Correlations introduced 
by such interactions can be understood in terms of ex- 
changeable processes. 

For dealing with such binary correlated systems the 
following representation was suggested by Laplace in 1774, 
and later rediscovered by de Finetti [23U5S] ; 



= / dyy n (l-y) 



N- 



i(y) 



dyg{y) = l . (5) 



This representation ensures Leibniz triangle rule and nor- 
malization, Yln=o ( 1 D r n = 1- Note that the non-negative 
g introduces the correlations in the stochastic process. In 
the case of independent variables, we have g(y) = S(y—p), 



i.e. r 



N 



N-n 



is the usual binomial factor for 
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independent trials. The fact that any stochastic process 
which is exchangeable to all scales N can be represented 
in this way is called the Laplace-de Finetti theorem [25l 
|27]|2"6"] (if processes are exchangeable only up to a scales 
Nmax positivity of g is not ensured and negative correla- 
tions can occur [26]). 

In order to pass from probabilities to distributions and 
from sums to integrals in a Riemann-Stieltjes sense in the 
thermodynamic limit N — > oo we use the well known prop- 
erties of the Beta function, B(a, b) = L dx x a ~ 1 (l — x) b_1 , 
to re-express the binomial factor in Eq. ([4]). Defining p^ = 
(N + l)p% turns Eq. g]) into 



^ dxx n {\- x) N - n g{x) 
Jq dxx n (l — x) 



N-n 



N 



^ N ■ 

n=0 



1 



As a consequence of this construction the natural domain 
of the limit distribution p(y) is the interval y G [0,1]. 
Using y% = (l/2 + n)/(iV+l), for n = 0, 1, . . . ,N, as dis- 
cretization of [0, 1] in the limit one gets l/(N + 1) — > dy 
and y„ — > y. Denote the natural number closest to the 
value (N + l)y - 1/2 by [(N + l)y] then the limit dis- 
tribution gets p(y) = limjv^oo p^ N+1 ^ y y These arguments 
imediately lead to 



p(y) 



lim 

N— >oo 



Jq dx [ X y{l-xf-y} N g{x) 
f'dx [ X v(l-x)i-v] N 



(7) 



Notice that the maximum of x v (l — x) 1 ^ with respect to 
x is obtained at x — y and gets amplified above all other 
values of the function by the power N. For N — > oo we 
therefore get a delta sequence and 



5(x — y) = lim 

N~*oc 



[x y {l~x) 



N 



J dx [xv(l - xy-y 



Consequently, 



p(y)= dx6{x -y)g(x) = g(y) 



(8) 



(9) 



i.e. limjv— ►<» 53™ \Pn ~ d{Vn)\ = 0- The limit distribution 
p and the function g generating the Leibniz-triangle 
are in fact identical! Therefore, once a desired limiting 
distribution on [0, 1] is given, one can simply write down 
the sequence of r„ which is generating it! Though it is 
known that g may become negative [26] , if scale- invariance 
is broken and the model is exchangeable only up to a max- 
imal scale N max , this interpretation of g makes it clear 
why non- negativity of g is required in the limit N — ► oo. 
The limit distribution p(y), i.e. g(y), of binary exchange- 
able processes is defined on y S [0,1], where y is ba- 
sically the ratio of events n/N in the limit of large N. 
However, prototypical processes, e.g. spin-systems or ran- 
dom walks, depend on binary processes (spins up/down 
or number of left/right steps). Yet, their associated ob- 
servable variable is not the ratio n/N of binary events 
but some other descriptive variable z (the magnetization 



or the position of a random walk). Since the domain of 
the descriptive variable z will in general not be [0, 1] (e.g. 
typically [—1,1] for spin systems and [—00,00] for random 
walks) it has to be clear from the beginning that the form 
of the limit distribution of the descriptive variable z will 
not only depend on g but also on the relation between 
z and the binary process, i.e. on y = n/N, which may 
or may not depend on the scale N. E.g., if equidistant 
spacing on intervals [— Rn,Rn] C [—00,00] as Rn — * 00 
is considered (as is the case of CLT) this relation will 
in general be a variable transformation depending on the 
scale N. For simplicity we only consider the case of N in- 
dependent transformations of variables between y and z. 
We assume a symmetric distribution G(z) of the descrip- 
tive variable z that can be obtained by a transformation 
of variables defined by dzG(z) — dyp(y). We will call a 
strictly monotonous increasing antisymmetric functions /, 
such that the effective stochastic variable z will take the 
values z„ = /(2j/„ — 1), a symmetric representation of 
the binary process. To be clear, z = f(2y — 1) is precisely 
the change of variables relating the generating function g 
of the binary process with the limit distribution G of the 
observable process z. In particular, G and g are one-to-one 
related by 2/'(2y- l)G(/(2y- 1)) = g(y). Each pair (/, g) 
exactly defines the distribution G of the binary process. 
Moreover, fixing G and the representation / uniquely de- 
termines g. Inversely, for a given observable distribution 
G, any pair (/, g) that represents G can serve as a stochas- 
tic model for G. A physical view on the meaning of / and 
g can be given by Galton's board, where / corresponds 
to the positions of needles on the board while g fixes the 
probabilities of balls being reflected to the left or right at 
the condition of hitting some needle. 

We now proceed to derive a stochastic model with 
q-Gaussian limit distributions. For the case of q < 1, 
the g-Gaussian G q (z) = -J-[l - (1 - q)z 2 } 1 / 1 ^ is de- 



fined on a compact support, \z\ < 



To identify it 



with the limiting distribution we have to map z to the 
unit interval by z — f(2y — 1) — ^==(2y — 1). Since 

the support of the q-Gaussian G q (z) for q < 1 is com- 
pact an affine transformation is the natural choice for 
the change of variables. Under an affine variable trans- 
formation G q (z) — * p(y) — -^-4y~ (1 — y) 1 ^ , where 

l) . Consequently, by introduc- 

1 

1, we get g(y) — tMikzli] 1 



Z q = 4B 



^- + 1 - L - 

l-q ^ X ' I-5 



ing the notation v 
and, finally, by using Eq. |5]), 



1 



B{v,v) 



B(u + n,v + N -n) 
B(v, v) 



(10) 



To retrieve the original g-Gaussian G q (z) one has to per- 
form the inverse coordinate transformation. This maps the 
discretization 



>£ = f(2y»-l) ) i.e., ^ = (2^-l)/0^(ll) 



and the discretization-width becomes dz = 2/(y/l — q(N- 
1)), which takes into account the the factor 2/'(2y^ r - 
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1) = 2/^/1 — q from the change of variables. Analogous 
to p the discrete formulation of the g-Gaussian G q reads 
(dz)-^)rf and 



F. 



N 



km "£\F» -G q (z»)\^0 



(12) 



This is exactly the model which was heuristically found in 
[24] , There the model was presented as 



A {N + v- 1)(N + v-2)...v 



(JV + 2i/-l)(iV + 2i/-2)...(2i/) ' 



(13) 



for integer values of v, which is obviously a particular 
case of Eq. (fTO]) . 

Since the g-Gaussian has no compact support for q > 1 
the situation becomes more involved, since now one has 
to map the real axis to the unit interval [0, 1]. Such a q- 
Gaussian G q (z) with z € [—00,00] may be thought of as 
distributions of the distance z a peculiar random walk has 
covered in the long time limit. There is no map / now that 
is a natural choice, a priori. However, in order to explicitly 
compute the probabilities one may choose a map such 
that is given in terms of Beta functions, as before. 
This leads to the situation that pairs of g-Gaussians, one 
with q < 1 and another with q > 1, can be generated by 
the same binary process and differ only in terms of the 
representation /. 

Let us find a map / such that f(2y — 1) again maps y € 
[0, 1] to z S [—00,00]. Using the normalization condition 
of the g-Gaussian this variable transformation implies 1 = 
f^dzGqiz) = J_ 1 dxf(x)G q (f(x)) and identify 

g(y) = 2f'(2y-l)G q (f(2y-l)) . (14) 

A particular / to compute in closed form is given by 



/(») = 



y 



1 



y/i - y 2 Vq 



1 



and f(y) 



(1 - v 2 ) 



-3/2 



(15) 



Noting, that this model / for the g-Gaussian implies [1 — 

(1 - q)f(y) 2 ] 1/{1 ~ q) = (1 - y 2 ) 1/(9_1) . Inserting this into 
Eq. (fT4|) we finally get the stochastic model, 



(?) 



B(p + n,p + N 
B(n,n) 



(16) 



where we have used the notation n = —^-r — i. All that is 

' q— 1 2 

left is to place everything correctly. We use the same dis- 
cretization y% = (1/2 + n)/(N + 1) for the interval [0, 1] 
as for q < 1. Again, dy„ = l/(N + 1) is the width of 
the discretization on [0, 1] and with p 1 ^ — (dy^) -1 ( n )r^, 
we get that limAr— «> £) n \p% - g(y%)\ = 0, as a power of 
TV as in the < q < 1 case. To retrieve the g-Gaussian 
G q (z), one has to perform the inverse coordinate trans- 
formation y £ [0, 1] — > z € [—00, 00]. This is a little more 
complicated since the discretization width now depends on 
N and n. In particular, mapping back to [—00, 00] gives 
the discretization z„ = f(2y% — 1), i.e. z^ 1 = (y^ — 



l/2)/y/{q - l)y%(l-y%), and the discretization width 
becomes dz% = [y% (1- y%))- 3/2 /(4{N + l)y/q~~l). The 
discretized version of the g-Gaussian for q > 1 now reads, 
F» = (dz^)- 1 fflr". In the limit we get lim^ £„ \F»- 
Gq{zn)\ — > 0, as a power of N. Comparing Eq. (jTO]) and 
Eq. (fl"T)|) it is obvious that for the two models of the q- 
Gaussian distribution (q < 1 and q > 1) the generat- 
ing binary processes are identical whenever i/ = /i. 
v and /1 are functions of q and g, and v(q) = p(q) es- 
tablishes a relation between two g-Gaussian distributions 
generated by an identical exchangeable binary stochastic 
7 ~ oq ■ q increasing from —00 to 1 yields 



process, i.e. q 



5-3q' 



q decreasing from 5/3 to 1. Therefore, the models of q- 
Gaussian distributions with q < 1 are conjugate with the 
model of g-Gaussian with q € [1,5/3] in the sense of be- 
ing driven by an identical binary stochastic process. The 
class of g-Gaussian distributions with g £ (5/3,3], which 
is exactly the class of normalizable g-Gaussian distribu- 
tions with diverging second moment, are not identified 
with a q < 1. The corresponding binary processes are 
unique in this sense. For instance, choosing 1 < g = 2 
requires 1 > g = 3, which is impossible. For q = 2 it 
follows that p = 1/2 and the associated binary process 
r N = b (1/2 + 71, 1/2 + N - n) /B (1/2, 1/2) has no rep- 
resentation for any g < 1, and r§ — ^g^ 2 /^ 2 /^' 1 — 
2-2JV ^2 \ course one can choose many families (f q , g q ) 
of models for g-Gaussian limit distributions and, for each 
strictly monotonous function g : [— 00, 1] 1— > [1,3], it is 
possible to construct conjugate families of models in the 
sense that they are generated by the same family of bi- 
nary processes. In fact when the function L q is defined by 

J^ q(z) dz'G q (z') = f* dz"G q{q) (z") and (/„ g q ) is a model 
for a g-Gaussian with g < then (L q (f q ), g q ) is a conju- 
gate model of the g-Gaussian with g = g(g). This allows to 
conjugate families of models for different dualities recog- 
nized in g-statistics, e.g. q = (5 — 3g)/(3 — g) for g-Fourier 
transforms [28] , 

We have explicitly derived two possible stochastic mod- 
els for correlated and exchangeable binary random vari- 
ables, which lead to exact g-Gaussians as the limiting dis- 
tributions, while strictly satisfying Leibniz rule, 



B 



2-g 
1-9 



n, 



2-9 

l-q 



N 



1/2 N 



B 



3-g 
2q-2 



( 2-9 2-g 
\l-q> l-q 



3-q 
2q-2 



N 



B 



( 3-g 3-g \ 
\^2g-2 ' 2g-2 J 



if q < 1 . 



if q = l. 



if 1< q < 3 , 



2g-2' 2g-2 

(17) 

5=1 can be obtained by both q — > 1^. Eq. (jTTJ) follows 
for the affine representation for g < 1 and for the represen- 
tation / given in Eq. (|15p for q > 1. For different families 
of representations f q of the q-Gaussians the equation has 
to be adapted accordingly. 

Lately the question whether Sq can be the extensive 
entropy for specific systems for Q < 1 has been discussed, 
e.g. in [9|fT0] . With the Laplace-de Finetti representation 
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of the observable variable of model is embedded on the 
real axis on each scale N. We have demonstrated this 
with the q-Gaussian distributions as an example. We have 
constructed scale invariant models for q-Gaussians for all 
q < 3 and have shown how notions of conjugate families 
of processes can be constructed. We have shown that the 
Boltzmann Gibbs entropy is the extensive entropy for all 
stochastic processes generated by scale-invariant binary 
processes. 

We acknowledge useful remarks by A. Rodriguez, V. 
Schwammle, E.P. Borges and E.M.F. Curado. R.H. and 
S.T. acknowledge great hospitality at CBPF and partial 
support by Faperj and CNPq (Brazilian agencies) and 
FWF Project P19132. 



Fig. 1. Nonadditive entropy Sq vs. N for a g-Gaussian with 
q = 1.5, for Q = 0.92, 1 and 1.02. Clearly, only Q = 1 is 
extensive (i.e., Si(N) oc N, for N > 1). 



we consider the entropy Sq for large N and prove that 
the value of Q for which Sq is extensive with respect to 
the distribution of binary states is Q = 1, i.e. q en t = 1, by 
straight foreward calculation: 



s Q [g] = 



i-EtoO^F 



((^v+i)(~)) 1 " Q 

Z^n=0 N+l 



J Q 1 dxx"(l-x) iv -"g(x) 



f£dxx"(l-x)"-r 



Sidy \s l dx\xy{i-xY-v] N 



1-0 

The first line follows from using the exchangeability-property 
together with the definition of the q-entropy. The second 
line can be obtained by inserting Eq. for and us- 
ing the properties of the /3-function. The third line follows 

for large AT, where Yln=o 1/(^+1) — * Jo 1 anc ^ usm S 
Eq. ([8|). To find the index Q where Sq gets extensive re- 
quires dS f^ to be a positive constant for N ^> 1 for some 
specific value of Q which is called g on t. We find, by using 
Stirling's approximation, 



dS Q [g](N) 



dN 



2 JV(i-Q) 5 (i/ 2 )hi2 



Hence, for any g(l/2) > 0, Sq(N) growing linearly with 
large N requires that Q = 1, since otherwise Sq(N) grows 
exponentially with N, i.e., q en t = 1 for all such models. It 
should be noted that particular choices of representations 
/ and the observable process (distribution functions) only 
result in particular choices of g. Since the result is valid for 
all g the result is basically independent of the particular 
choice of representing an exchangeable process. 

We have shown that the Laplace-de Finetti represen- 
tation is a suitable framework which allows to generate 
scale invariant probabilistic models by fixing their limit 
distribution and a map which describes how the domain 



References 

1. Anteneodo C. and Tsallis C, J. Math. Phys., 44 (2003) 
5194 . 

2. Fuentes M.A. and Caceres M.O., Phys. Lett. A, 372 (2008) 
1236-1239. 

3. Borland L., Phys. Lett. A, 245 (1998) 67. 

4. Plastino A.R. and Plastino A., Physica A, 222 (1995) 347. 

5. Tsallis C. and Bukman D.J., Phys. Rev. E, 54 (1996) 
R2197. 

6. Umarov S., Tsallis C. and Steinberg S., Milan J. Math., 76 
(2008) [DOI 10.1007/s00032-008-0087-y]; for a more ped- 
agogical version, see Queiros S.M.D. and Tsallis C, AIP 
965, 21 (New York, 2007). 

7. Tsallis C, J. Stat. Phys., 52 (1988) 479; Introduction to 
Nonextensive Statistical Mechanics - Approaching a Com- 
plex World (Springer, New York, 2009), in press. 

(18) 

8. Penrose O., Foundations of Statistical Mechanics: A De- 
ductive Treatment (Pergamon, Oxford, 1970), p. 167. 
Tsallis C, Gell-Mann M. and Sato Y., Proc. Nat. Acad. 
Sci., 102 (2005) 15377. 

10. Caruso F. and Tsallis C, Phys. Rev. E, 78 (2008) 021101. 

11. Pluchino A., Rapisarda A. and Tsallis C, Eur. Phys. Lett., 
80 (2007) 26002; Physica A, 387 (2008) 3121. 

12. Anteneodo C. and Tsallis C, Phys. Rev. Lett., 80 (1998) 
5313. 

13. Yamaguchi Y., Bouchet G. and Dauxois T., JSTAT, (2007) 
P01020. 

14. Bouchet G. and T. Dauxois T., Phys. Rev. E, 72 (2005) 
045103 (R). 

15. Upadhyaya A., Rieu J.-P., Glazier J. A. and Sawada Y., 
Physica A, 293 (2001) 549. 

16. Thurner S., Wick N., Hanel R., Sedivy R. and Huber L.A., 
Physica A, 320C (2003) 475-484. 

17. Daniels K.E., Beck C. and Bodenschatz E., Physica D, 193 
(2004) 208. 

18. Burlaga L.F. and Vinas A.F., Physica A, 356 (2005) 375. 

19. Burlaga L.F., Vinas A.F., Ness N.F. and Acuna M.H., As- 
trophys. J., 644 (2006) L83-L86. 

20. Douglas P., Bergamini S. and Renzoni F., Phys. Rev. Lett., 
96 (2006) 110601. 

21. Liu B. and Goree J., Phys. Rev. Lett., 100 (2008) 055003. 

22. Tirnakli U., Beck C. and Tsallis C, Phys. Rev. E, 75 
(2007) 040106(R); 0802.1138 [cond-mat.stat-mech]. 

23. Grassberger P., arXiv:0809.1406 [cond-mat.stat-mech]. 



(19) 



6 



Rudolf Hanel et al.: Title Suppressed Due to Excessive Length 



24. Rodriguez A., Schwammle V. and Tsallis C, JSTAT (2008) 
P09006. 

25. de Finetti B., Atti della R. Academia Nazionale dei Lin- 
cei 6, Memorie, Classe di Scienze Fisiche, Mathematice e 
Naturale, 4 (1931) 251299; Heath D., Sudderth W., Am. 
Stat., 30 (1976) 188-189. 

26. Jaynes E.T., in Bayesian Inference and Decision Tech- 
niques with Applications, Goel P.K. and Zellner A., eds. 
(North-Holland, Amsterdam, 1986), p. 31. 

27. Hewitt E., Savage L.J., Trans. Am. Math. Soc., 80 (1955) 
470-501. 

28. Nelson K.P.. Umarov S.. llarXiv:0811.3777H . 



